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ABSTRACT 

Aims. The effects of multi-layered clouds in the atmospheres of Earth-like planets orbiting different types of stars are studied. The 
radiative effects of cloud particles are directly correlated with their wavelength-dependent optical properties. Therefore the incident 
stellar spectra may play an important role for the climatic effect of clouds. We discuss the influence of clouds with mean properties 
measured in the Earth's atmosphere on the surface temperatures and Bond albedos of Earth-like planets orbiting different types of 
main sequence dwarf stars. The influence of clouds on the position of the habitable zone around these central star types is discussed. 
Methods. A parametric cloud model has been developed based on observations in the Earth's atmosphere. The corresponding optical 
properties of the cloud particles are calculated with the Mie theory accounting for shape effects of ice particles by the equivalent 
sphere method. The parametric cloud model is linked with a one-dimensional radiative-convective climate model to study the effect 
of clouds on the surface temperature and the Bond albedo of Earth-like planets in dependence of the type of central star. 
Results. The albedo effect of the low-level clouds depends only weakly on the incident stellar spectra because the optical properties 
remain almost constant in the wavelength range of the maximum of the incident stellar radiation. The greenhouse effect of the 
high-level clouds on the other hand depends on the temperature of the lower atmosphere, which is itself an indirect consequence 
of the different types of central stars. In general the planetary Bond albedo increases with the cloud cover of either cloud type. An 
anomaly was found for the K and M-type stars however, resulting in a decreasing Bond albedo with increasing cloud cover for certain 
atmospheric conditions. Depending on the cloud properties, the position of the habitable zone can be located either farther from or 
closer to the central star. As a rule, low-level water clouds lead to a decrease of distance because of their albedo effect, while the 
high-level ice clouds lead to an increase in distance. The maximum variations are about 15% decrease and 35% increase in distance 
compared to the clear sky case for the same mean Earth surface conditions in each case. 
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1. Introduction 

Cloud particles can have an important impact on the climate of 
planetary atmospheres by either scattering the incident stellar ra- 
diation back to space (albedo effect) or by trapping the infrared 
radiation in the atmosphere (greenhouse effect). The answer to 
the question which of these effects dominates for a given cloud 
type depends on a variety of cloud parameters, the most impor- 
tant of which are the cloud composition (water, carbon dioxide 
etc.), the size distribution of the cloud particles, the optical depth 
of the cloud layer, multi-layered cloud coverage and the cloud 
altitude. 

In the case of the well known Earth atmosphere, where 
clouds are a very common phenomenon with a mean global 
cloud coverage of more than 50%, low-level water clouds have a 
net cooling effect on the surface, while high-level ice clouds ex- 
hibit a greenhouse effect, resulting in surface heating. A compre- 
hensive review about the climatic effe cts of c l ouds i n the Earth 
atmosphere can be found in e.g. iKondratvevi d!999l) and refer- 
ences therein. The albedo and the greenhouse effect are directly 
correlated with the wavelength-dependent optical properties of 
the cloud particles. The incident stellar spectra in combination 



Send offprint requests to: D. Kitzmann 



with these optical properties of the cloud may therefore play an 
important role for the surface temperature. 

Model calculations regarding the habitability of planets or- 
biting dif ferent types o f cent ra l stars have a l ready been done 
bv e.g. iKasfing et alJ (Il993h . ISegura et alJ d2003l I2005I) or 
Selsis et alJ (I2007I) . Still, these models lack a detailed treatment 
of the cloud radiative forcing in the radiative transfer and aim 
only to mimic the cooling effect of low-level clouds by an ad- 
justed surface albedo. Focuss ing on the inn er boundary of the 
habitable zone around the Sun Kasting ( 1988) included the effect 
of one very extended water droplet cloud with rather simplified 
optical properties in some of the model scenarios. Without tak- 
ing clouds explicitl y into account in t heir a tmospheric climate 
model calculations Kalt enegger et alJ (12007b studied the emis- 
sion spectra of Earth-like planets at different evolutionary stages. 

We study the effects of different incident stellar spectra in 
conjunction with multi-layered clouds of different types on the 
surface temperatures of Earth-like extrasolar planets, which de- 
termine the potential habitability of a terrestrial planet. As a 
first approximation, a simplified cloud description scheme us- 
ing parametrised size distribution functions is used. Other cloud 
properties like optical depth and cloud top pressure have been 
taken from measurements. The parametric cloud model is de- 
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Table 1. Mean cloud properties obtained from surface and satellite measurements. 



Cloud type 


Cloud coverage (%) a 
Ocean Land Global mean 


Global average cloud properties b 
Top temperature Top pressure Optical depth c 


Low-level clouds 


47 


22 


39.5 


281.1 K 0.826 bar 4.7 


Stratus 


12 


5 


9.9 




Stratocumulus 


22 


12 


19.0 




Cumulus 


13 


5 


10.6 




Cumulonimbus 1 ' 


6 


4 


5.4 




High-level clouds 


12 


22 


15.0 


227.5 K 0.267 bar 2.2 


Cirriform 


12 


22 


15.0 





a taken from surface observations (Warren et al. 2007), the global mean coverages have been derived by a weighted average from the coverages 

over ocean (70%) and land (30%) 

b from satellite observations within the ISCCP (see lRossow & Schifferl {[999) for details) 
c optical depth at 0.6 yum 

d not included in the model and in the calculation of the mean coverages 



scribed in Sect. [2] To study the basic climatic effects, our cloud 
scheme has been coupled with a one-dimensional radiative- 
convective climate model which includes the possibility to ac- 
count for different amounts of cloud coverages as well as the 
partial overlap of two cloud layers, as described in Sect. [3] In 
Sect. [4] we apply this climate model including the parametrised 
cloud description to the modern Earth atmosphere to verify the 
applicability of our model approach. Resulting surface temper- 
atures and Bond albedos for different cloud types in the atmo- 
spheres of Earth-like planets orbiting different types of central 
stars and implications for the positions of the habitable zones 
are presented in Sect. [5] 

2. Cloud model description 

Based on the extremely well-studied properties of different 
cloud types occurring in the Earth's atmosphere we developed 
a (parametrised) multi-layered cloud description scheme. Two 
different kinds of cloud layers are considered here: low-level 
water droplet and high-level water ice clouds. The correspond- 
ing global and temporal average cloud coverages resulting from 
long-term su r face observations have been published by e.g. 
IWarren et al.1 (120071) and are given in Table Q] The measured 
cloud coverages already include a certain amount of overlap be- 
tween the different c loud l ayers, which was not further speci- 
fied bv l Warren et al.l (120071) . Average properties of Earth clouds, 
for instance the cloud top temperature and pressure as well as 
their optical depth have been derived by long-term satellite- 
based measurements within the International Satellite Cloud 
Climatology Project (ISCCP) for example. Global and temporal 
mean clo ud properties based on th ese surveys have been pub- 
lished bv lRossow & Schifferl (1 19991) and are also summarised in 
Table □ 

Other cloud types present are not taken into account. With 
a mean coverage of 17% the most common mid-level clouds, 
altocumulus clouds, have been reported to be radiatively neu- 
tral, which mean s that their albedo and green house effect bal- 
ance each other dPoetzs ch-Heffte ret al.l Il995h . which justifies 
our approach to neglect them. Cumulonimbus clouds are also 
excluded from our cloud description scheme. These clouds can 
extend up to 10 km, which makes it difficult to include them in 
our one-dimensional climate model (see Sect. [3]). For such ex- 
tensive vertical clouds not only the radiation entering the cloud 
from the top or bottom has to be considered, but also photons 



entering the cloud sideways have to be taken into account. This 
though is essentially a three-dimensional problem which cannot 
be treated by a one-dimensional model. Since such clouds can 
not be properly included and also present only 5% of the global 
cloud coverage, we decided against implementing them into our 
column model. 

Apart from these measured properties the size distribu- 
tion of the cloud particles has also to be known to calculate 
the wavelength-dependent absorption and scattering properties 
of a cloud. Usually it is not appropriate to descr i be a c loud 
with uniform-size particles (but see e.g. iMitchelll (120021) and 
McFarquhar & Heymsfield (1998) concerning the applicability 
of an effective radius). The radius a of a cloud particle has there- 
fore to be treated as a random variable, characterised by a dis- 
tribution function f(a). Normally these distributions have to be 
calculated from first principles ac counting for all relevant m i- 
crophysical cloud processes (see Prupp acher & Klettl (Q~997)), 
which involves a great deal of computation. For the Earth's at- 
mosphere it is possible though to derive parametric analytical 
distribution functions based on measurements. Since the focus 
of this work is on Earth-like planetary atmospheres, we assume 
that these analytical distribution functions are also valid for the 
atmospheres considered here. This assumption neglects any pos- 
sible influence of different atmospheric properties (like differ- 
ences in the atmospheric dynamics or chemical composition) on 
the distribution function of the cloud particles. Such differences 
could arise when considering planets with for instance different 
rotation periods or a different landmass distribution compared to 
Earth. Therefore the cloud parametrisations used here are limited 
to Earth-like planets. 

2. 1 . Low-level cloud particle size distribution 

Observations of low-level clouds in the Earth's atmosphere show 
that the measured size distribution of the cloud particles can be 
well-represented by a log-normal distribution 

tt \ 1 / Qna-lna m ) 2 \ 

Measured parameters for the particle density n, the mean particle 
ra dius g m and the standa rd deviation <x can be found for instance 
in Kokhanovskj] (2006). The mean values for continental and 
maritime water clouds are given in Table [2] 



D. Kitzmann et al.: Clouds in the atmospheres of extrasolar planets. I. 



3 



Table 2. Parameter sets for the description of maritime and con- 
tinental water clouds by a log-normal distribution. 



Parameter 


Maritime clouds 


Continental clouds 


n (cm -3 ) 


91.0 


254.0 


a m (K) 


6.0 


4.0 


cr 


0.4 


0.4 



The distribution functions using these sets of parameters are 
averaged according to the total ocean (70%) and land (30%) sur- 
face area of the Earth. The resulting mean distribution function 
is used to represent low-level clouds in our climate model. We 
also assume that all low-level cloud droplets are spherical parti- 
cles composed of pure liquid water. The optical properties from 
the obtained distribution of the cloud particles can then be cal- 
culated with the Mie theory (see Sect. I2.31 >. 



2.2. Size distributions of high-level cloud particles 

High-level clouds represent a much more complex system, be- 
cause the ice crystals can have a huge variety of shapes, which 
makes the description by a single distribution function and the 
calculation of the optical properties much more complicated. 
The most common shapes present in these clouds are solid and 
hollow columns, plates, bullets and bullet rosettes. For simplic- 
ity all ice particles are considered as solid hexagonal columns 
throughout this work. 

Based on in - situ measurements in cirrus clouds, 
iHevmsfield & Plattl d 19841) derived analytical size distributions 
for high-level ice clouds using a power law size distribution 



f(a c ) = Aa B c 



(2) 



where A is the intercept parameter, B the slope, and a c the col- 
umn's maximum dimension. Depending on the cloud's temper- 
ature, a bimodal or unimodal distribution must be used. For the 
temperature range in our case (see Table [T]) a unimodal distribu- 
tion is sufficient to describe the measured size distribution. The 
parameters A and B are temperature dependent and for the tem - 
perature range in question given by (Hevmsfield & Piatt 1984) 



2.55 ■ 10- 

A = 5— 

100 B 

B = -3.23 . 



//m 



(3) 
(4) 



Sin ce the size distribution derived by IHevmsfield & Plattl 
(1 19841) depends only on the maximum particle dimension, we 
additionally need to prescrib e the aspect ratio of the particles to 
fully describe the columns. Hevmsfieldi ( 1 1 9721) provided analyti- 
cal expressions for the aspect ratio T of the column base width D 
and the column length L valid for temperatures near m -20 °C. 
For L < 200 fim the aspect ratio is given by 



T(L) = 2 
whereas 
T(L) = 5 ■ L 1 



(5) 



0.59 



is used for crystals with L > 200 /im . 

Because the length L is always longer than the base width 
D, for columns with these aspect ratios the maximum dimension 
a c , which is used to parametrise the distribution function (Eq.|2]i, 
refers to the crystal's length in all cases. 



Even assuming exclusively solid hexagonal columns as par- 
ticle shapes, the calculation of the corresponding optical prop- 
erties is nevertheless more complex than in the case of spheri- 
cal droplets because the Mie theory cannot be applied directly. 
There do exist some applications for the derivation of the opti- 
cal properties of these non-spherical particles, but they are either 
limited in their validity range (e.g. geometrical optics or ray trac- 
ing methods) or need an excessive amount of computing time 
(e.g. finite-time-domain theory). 

However, it is possible to introduce so-called equivalent 
spheres, i.e. the non-spheric particles are replaced by spheres, 
which mimic their optical properties. In this way the Mie the- 
ory can be used again allowing a fast computation. Commonly 
used are equal-volume spheres and equal-area spheres, whereby 
a sphere has the same volume or the same surface area as the 
respective non-spherical particle. In both cases the number den- 
sity of the equivalent spheres and the non-spherical particles 
remain the same, while either conserving the total volume (in 
the equal-volume sphere case) or the total surface area (for 
the equal-area sph ere case) of the non-spher ical particles. Still, 
as pointed out by iGrenfell & Warren] d!999l) . these approaches 
yield mostly too small scattering albedos and too large asymme- 
try parameters in comparison with the non-spherical particles. 
A much better agreement is achieved by using spheres having 
the same volume-to-surf ace ratio as the non-spherical particles 
( IGrenfell & Warrenl 1999). But to conserve the total volume and 
area of the non-spherical particles, the number density of the 
equivalent spheres has to be adapted. Compared to the equal- 
volume or area spheres, the sizes resulting from the volume-to- 
surface equivalent sphere method are generally smaller, which 
in turn leads to smaller asymmetry parameters and larger scatter- 
ing albedos. This implies that the volume-to-surface equivalent 
sphere approach offers a much better approximation for the cal- 
culation of the optical properties of the non-spherical particles 
in most cases. 

The application of the volume-to-su rface sphere approac h 
for hexagonal columns was published by Neshvba et al. (2003). 
Following their treatment the radius of the equivalent spheres for 
hexagonal columns is given by the expression 



3rV3 



(7) 



(8) 



4r+ V3 

while the number density of the spheres n s is determined by 
V3 

As noted by Neshvba et al. (2003) these approximations 
have only been tested for use in energy budget studies, i.e. when 
using angle-averaged properties of the radiation field (e.g. the 
spectral flux) as done in this study (cf. 13.21 1. Their applicabil- 
ity might though be limited concerning the calculation of angle- 
dependent spectral intensities. In particular the required assump- 
tion of randomly oriented columns is questionable, because most 
of the columns are falling with their long axes paralle l to the 
ground as confirmed by observations of e.g. lOnol d!969l) . 



(6) 2.3. Optical properties of cloud particles 



For spherical particles within the size interval (a + da) the usual 
transport coefficients are given by 



Xx{d)da = na 2 Q e ff(a)da 
~s,i(a)da - no 2 f (a)da 



(9) 
(10) 
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ncP'QY' f(a)da - XA(fl)da - sx(a)da . 



(11) 



where \a denotes the extinction coefficient, the (total) scat- 
tering coefficient and the absorption coefficient. Integration 
over the whole size range yields the complete contribution from 
the distribution function f(a) to the total values of the optical 
properties 



Xa 



-r 

Jo 

-f 

Jo 



X,\(a)da 



sx{a)da 



*vt = I Kx(a)da . 
Jo 



(12) 
(13) 
(14) 



The required optical properties (extinction efficiency Q e , 



scattering efficiency (7j ca , and absorptio n efficiency Qjf s ) are cal - 
culated here using the Mie theory dBohren & Huffmanlll983h . 
Since the Mie calculations become computationally prohibitive 
for large size parameters, the large particle limit is employed 
in these cases. Using the large particle limit the efficiencies are 
given by the simple expressions 



Q 



= 2 
= 1 + 



Qf 



i - Qf 



(15) 
(16) 
(17) 



where the reflection efficiency Q r ^ n can be calculated with the 
geometric optics approximation which is valid for particles 
much larger than the considered wavelength. 

Instead of the full scattering phase function p,\(6,a), the 
asymmetry parameter ^(a), which is defined by 



gA(a) 



1 r 

2 Jo 



Pa(0, a) cos sin 9d6 



(18) 



for a single particle of radius a and the scattering angle 6 can be 
used to approximate the effects of multiple scattering. The asym- 
metry parameter can be written as a function of the Mie c oeffi- 
cients and the size parameter (see lBohren & H uffman 19831) . 

For a continuous distribution function as in the cloud de- 
scription, the average value 



8a 



na 2 gA (a)QTf(a)da 
f Q °° na 2 Q s ™f(a)da 



(19) 



is used, applicable for multiple scattering. 

In the case of the low-level clouds we use the refractiv e 
indices for pure liquid water taken from ISegelsteinl (Il98ll) . 
whereas f or the high-level c l ouds t he indices for water ice, pub- 
lished by Warren & Brandt (2008) are adopted. The refractive 
indices are assumed to be independent of temperature. 

The results of the Mie theory calculations for high- and 
low-level cloud size distributions have been scaled according 
to the measured optical depth from Table [1] and are shown in 
Fig. Q] Scattering dominates the radiative effects in the wave- 
length range below 1 //m, where the maximum of the incident 
stellar spectra is located. With a scattering albedo of nearly 
tox ~ 1, absorption is negligible in this wavelength range. Thus 
the clouds will mostly scatter the incident light at short wave- 
lengths. However, due to high asymmetry parameters (gj ss 0.8) 
much of the stellar radiation will be scattered in the forward di- 
rection, i.e. it will still reach the planetary surface, and only a 
small part is reflected back to space. In the longer wavelengths 




Fig. 1. Calculated optical properties of the cloud model. Upper 
diagram: optical depth of the low-level water cloud (solid line) 
with the individual contributions of absorption (dotted line) and 
scattering (dashed line); upper middle diagram: optical depth 
of the high-level ice cloud (solid line) with the individual contri- 
butions of absorption (dotted line) and scattering (dashed line); 
lower middle diagram: asymmetry parameter g,\ of the water 
cloud (solid line) and ice cloud (dashed line); lower diagram: 
scattering albedo co^ of the water cloud (solid line) and ice cloud 
(dashed line). 

range of the outgoing thermal radiation around 10 fim, the radia- 
tive effects are much more complex, because the optical prop- 
erties, especially the absorption and scattering optical depths, 
show large variations. The different effects for this case are ex- 
plained in detail in Sect lS. 1 1 

3. Radiative-convective climate model 

3. 1 . Basic assumptions 

For the atmospheric model calculations we use a one- 
dimensional radiative-convective climate model, which is based 
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Table 3. Properties of the different central stars. 



Fig. 2. Atmospheric chemical profiles of O3 (solid line), CH4 
(dashed line) and N2O (dotted line). The profiles have been de- 
rived from a radiative-convective climate model coupled with 
a photoche mistry model r eprese nting the modern Earth atmo- 
sphere (see lGrenfell et al.1 (120071) ). 



on th e model developed and described bv lKasting et al.l (fl984) 
and Pavlov et al. (2000). In their atmospheric model the impact 
of clouds is not explicitly treated. The effects of clouds are only 
taken into account indirectly by adjusting the value of the plan- 
etary surface albedo to mimic the influence of clouds in the tro- 
posphere. Using the developed cloud model description (Sect.[2jl 
we include the climatic effects of multi-layered clouds directly 
into the climate model to determine for instance the radiative 
feedbacks of clouds on the surface temperature. Chemical feed- 
backs of clouds are not included yet. The atmospheric profiles of 
the major c hemical species obtained with a detailed photochem- 
ical model dGrenfeil et al. 2007) representing the modern Earth 
atmosphere are used. Since N2, O2, and CO2 are well mixed 
within the atmosphere, their atmospheric profiles are given by 
isoprofiles with mixing ratios of 78%, 21% and 0.0355% re- 
spectively. The profiles of CH4, O3 and N2O have been derived 
from the photochemical model for the modern Earth. This atmo- 
spheric composition is assumed for all calculations, thereby ne- 
glecting any change of processes influencing the chemical com- 
position of the planetary atmospheres, like different CO2 levels 
due to changes in silicate weathering, for example. 

For the relative humidity in the troposphere, from which the 
water profi le is calculated, the empirica l relative humidity dis- 
tribution of M anabe & Wetheraldl d 1967b is used. Within the tro- 
posphere the temperature is assumed to follow a moist adiabate, 
otherwise the temperature is calculated by the condition of ra- 
diative equilibrium. The measured Earth surface albedo of 0.13 
is applied, taken from globally averaged satellite measurements 
of the ISCCP. The surface emission is treated as black-body ra- 
diation determined by the surface temperature. 

3.2. Radiative transfer 

The radiative transfer within the climate model is split into two 
wavelength regimes: the stellar part, dealing with the wavelength 
range of the incident stellar radiation and the infrared part for the 
treatment of the thermal radiation wavelength range. 

The radiative transfer in the stellar part consists of 38 spec- 
tral intervals between 0.238 /im and 4.55 /im. The plane-parallel 
equation of ra diative transfer i s solved by a <5-two-stream quadra- 
ture method dToon et al.ll 1 989h . Gaseous absorption of 63, O2, 
H2O, CO2, and CH4 is treated with four-term correlated-k coef- 



Name 


Type 


T ef f (K) 


d (pc) 


a (AU) 


cr Bootis a,b 


F2V 


6722 


15.5 


1.89 


Sun c 


G2V 


5777 





1.0 


e Eridani a d 


K2V 


5072 


3.2 


0.61 


AD Leo c f 


M4.5V 


3400 


4.9 


0.15 



Habi ng et al J d200lh : b ICenarro et al 
Santos etal.U2004h : e iLeggett et ail 



ficients (cf. lSegura et~al1(l2003l) ). For 2 , C0 2 and N 2 , Rayleigh 
scattering is also considered. 

For the IR part a hemispheric-mean two-stream method is 
used, since the 5-two-stream method used in the stellar part can- 
not be used for th ermal radiation due to numerical inaccuracies 
dToon et al.lll989l) . The IR radiative transfer uses 16 spectral in- 
tervals between 3 yum and 1000 /im. 

Gaseous absorption is treated with the correlated-k method 
incorporated within the rapid radi ative transfer model (RRTM) 
developed by Mla wer et al.l (1 19971) . Considered species in the 
RRTM for gaseous absorption in the IR wavelength range are 
O3, N 2 0, H 2 0, C0 2 , and CH 4 . Scattering is neglected for 
gaseous species at these wavelengths. The same approach has 
been used for instance by Segura et al.l d2003l) . 

Note that the k-distributions used in the RRTM are only tabu- 
lated over a limited pressure and temperature range. In particular 
the temperature range is limited to within + 30 K of an Ear t h mid- 
latitude summer temperature profile (see iMlawer et alj (1 19971) 
for details). Beyond the tabulated range extrapolation is used to 
derive the k-coefficients, which might cause inaccuracies in the 
upper parts of the atmospheric temperature profiles, especially 
for situa tions strongly deviating fr om (mean) Earth conditio ns 
(see also lSegura et all d2003l 120051) and lvon Paris et al.l d2008l) 1. 



3.3. Incident stellar spectra for typical M, G, K, and F-type 
stars 

Atmospheric calculations of Earth-like planets around different 
main sequence dwa rf stars, M, F, G, and K- type stars, been have 
publi shed by e.g. iKasting et ail d 19931) or ISegura et all (|2003, 
2005). In order to be comparable we investigated the same stars 
in this study as representatives for the different stellar types. The 
sample of stars used for the calculations are the F-dwarf cr Bootis 
(HD 128167) as a typical F-type star, the Sun as the G-type star, 
the young active K-type star e Eridani (HD 22049), and the M- 
type dwarf star AD Leo (GJ 388). The basic stellar parameters 
(stellar type, effective temperature T e ff, distance of the consid- 
ered star to the Sun d, and distances from the planets to their 
host stars a, for which the stellar flux matches the solar constant) 
with corresponding references are shown in Table [3] 

High resolution spectra for the M-dwarf and the F-type 
star were taken from lSegura et al] (I2003L |2005). According to 
Segu ra et~ai1 d2003l) the F-type star spectrum is a composi- 
tion from IUE satellite spectra of <x Bootis between 115 nm 
and 335 nm and a synthetic spectrum derived from the well- 
established stellar at mosphere model of Kurucz dKuruczl ll979: 
iBuser & Kuruczll 19921). More details on the F-star spectrum can 
be found in lSegura et al.l d2003l) . The M-type star high resolution 
spectrum is a composite of IUE satellite dat a between 115.1 nm 
and 33 4.9 nm and an optical spectrum from Pette rsen & Hawlevl 
(1989) between 335.5 nm and 900 nm. It was extended to 
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Fig. 3. Incident stellar spectra for different central stars. Each ra- 
diation flux is scaled to the distance where the total energy input 
at the top of the planetary atmosphere equals the solar constant. 

2409 nm in the near infrared using spectra by iLeggett et al.l 
(1996) and a synthetic photospheric spectrum from the stellar at- 
mosphere model NextGen beyond 2410 nm. For furth er details 
on the M-star spectrum we refer to lSegura et al.l d2005l) . 

The K-type star spectrum is composed of IUE satellite data 
from e Eridani between 115 nm and 355 nm and a synthetic 
NextGen spectrum, taken from the grid of stellar atmosphere 
models of France Allard (http://perso.ens-lyon.fr/france.all ard/). 
The hi gh resolution spectrum of the Sun is taken from Guevmard 
(2004). This updated compilation of stellar spectra has been used 
in this study. 

A total integrated solar radiation flux of 1366 WirT 2 has 
been derived by i ntegration of the high resolution spectrum of 
Guevmard (2004) from 0.5 nm to 10 00 /mi. In contrast to the 
approach of ISegura et al.l d2003l 12005). the distances of the four 
representative Earth-like extrasolar planets to their respective 
central stars have been determined in this work in a way that 
the integrated incident stellar flux equals this solar constant. The 
corresponding orbital distances of the planets are shown in Table 
[3] and the stellar spectra of all four central stars incident at the 
top of the planetary atmospheres are shown in Fig. [3] 

Furthermore, all four spectra have been binned to obtain the 
integrated radiation fluxes required for the 38 spectral intervals 
from 0.238 /mi to 4.55 /mi used in the radiative transfer (cf. Sect. 
E2J. 
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Fig. 4. Temperature-pressure profiles of the Earth reference 
models. The solid line represents the clear sky case, the dotted 
line the model with the mean Earth cloud cover. The horizontal 
dashed lines represent the position of both cloud layers. 



layer, partial overlapping clouds etc.) and the clear sky case (cs) 
the radiative transfer is solved separately. The mean radiative 
flux F A is then determined by averaging all separately calculated 
fluxes F x,i weighted with the corresponding cloud coverage x-, 

i i 
F A = (1 - x) ■ F AiCS + 2^ F AJ Xj with x = ^ *i < 1 • (20) 

4. Earth reference model 

In order to verify the applicability of our cloud description we 
first performed model calculations for the modern Earth atmo- 
sphere. ISCCP measurements report a global mean Earth surface 
temperature value of 288.4 K, which we considered as reference 
below. Calculations for the clear sky case and for the measured 
Earth mean cloud cover (39.5 % low-level and 15 % high-level 
cloud cover) have been done also considering a partial overlap 
of ss 7% between the two cloud layers. The calculated values for 
the surface temperatures, the temperatures at the positions of the 
cloud layers, and the Bond albedos are summarised in Table |4] 
and compared to measured values taken from ISCCP data. 



3.4. Cloud-climate scheme 

The climate model used in this work allows for a multi-layered 
cloud structure: two different cloud layers and the possibility of 
partial overlap between both layers. 

The optical properties of the clouds are calculated accord- 
ing to the description in Sect. 12.31 and the results summarised 
in Fig. Q] The optical properties are assumed to be the same 
for all model atmospheres, i.e. no feedback of the atmospheric 
properties (e.g. different temperature profiles) onto the clouds is 
considered here. The altitude of both cloud layers is iteratively 
adjusted to match the measured pressure values shown in Table 
ffl 

To account for the radiative effects of clouds, their optical 
properties (optical depths, asymmetry parameter, and scattering 
albedo) have been introduced into both parts of the radiative 
transfer scheme. In order to account for different amounts of 
coverages and their partial overlap of multi-layered clouds in our 
model we developed a flux-averaging procedure. For every dis- 
tinct cloud configuration i (e.g. a low-level or a high-level cloud 



4.1. Surface temperatures 

For the clear sky case, the calculated surface temperature is 
about 5 K too high. Using the global mean cloud coverages 
of 39.5% for the high-level cloud and 15% for the high-level 
cloud (see Sect.|2]and TableQ}, the resulting surface temperature 
(287.7 K) agrees much better with the observed value. As we 
pointed out in Sect. [2] these measured cloud coverages already 
include a certain amount of overlap between the cloud layers, 
but this value was not specified by the measurements. We deter- 
mined a partial overlap of about 7% required to reach a surface 
temperature of 288.4 K. 

The temperature pressure profile of the clear sky case and 
of the cloudy case (without overlap) are shown in Fig. [4] The 
positions of the two cloud layers are denoted by the two dashed 
lines. The profile of the model including partial overlap is not 
shown because it shows no noticeable difference compared to 
the non-overlapping case. The profiles again indicate that the 
clear sky calculation results in too high temperatures, while the 
cloudy case resembles mean Earth conditions. The presence of 
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Value 


Clear sky 


Earth mean cloud cover 


including partial overlap 


measured Earth values 


Surface temperature 


293.3 K 


287.7 K 


288.4 K 


288.4 K 


Low-level cloud temperature 


287.3 K 


280.7 K 


281.4 K 


281.1 K 


High-level cloud temperature 


229.2 K 


220.1 K 


220.8 K 


227.5 K 


Bond albedo 


0.15 


0.27 


0.26 


0.3 



clouds changes the temperature profiles up through the tropo- 
sphere, while the temperatures in the upper atmosphere are al- 
most unaffected as expected. 

The temperatures of the low-level clouds also agree very well 
with the measured values. In contrast to that, the high-level cloud 
temperatures show a dif ference of about 6.7 K in comparison to 
the values published by Ross ow & Schifferl(ll999l) (see Table[TJi. 
Still, as was pointed out in their work, the cloud temperatures 
were derived neglecting cloud IR scattering. This resulted in an 
overestimate of their high-level cloud temperatures by a few de- 
grees so that the actual deviations of our calculated cloud tem- 
peratures are much smaller than 7 K. Because the optical depth 
and the coverage of the low-level clouds is larger than that of 
high-level clouds, clouds in the present Earth atmosphere have a 
net cooling effect as confirmed by our model findings. 



4.2. Bond albedo 

The calculated value of 0.15 for the planetary Bond albedo (see 
Table 2]i is much too small in the clear sky case compared with 
the observed Earth value of about 0.3. With clouds included in 
the model, the resulting Bond albedo of 0.27 (0.26 in case of 
overlapping clouds) agrees much better with the observed value. 
The remaining small difference can be explained by our neglect 
of the mid-level and cumulonimbus clouds, which would also 
contribute to the Bond albedo. 



4.3. Radiative flux profiles 

For a better understanding of the cloud radiative forcing, the up- 
ward and downward radiation flux-pressure profiles in the stellar 
and IR wavelengths range are shown in Fig. [5] 

The downward stellar flux clearly indicates the albedo effect 
of both cloud layers. Due to its larger optical depth and bigger 
coverage, the albedo effect of the low-level cloud is much more 
pronounced than that of the high-level cloud. The upward in- 
frared flux shows the blocking of the thermal radiation by the 
different cloud layers, i.e. the resulting greenhouse effect. The 
low-level cloud has almost no effect on the outgoing IR radia- 
tion and therefore exhibits no noticeable greenhouse effect. This 
yields a net albedo effect (see the steps in the flux profiles in Fig. 
01, which leads to a cooling of the lower atmosphere and also 
of the planetary surface. The high-level cloud on the other hand 
traps more IR radiation in the lower atmosphere than it blocks 
the incident solar radiation. Thus, for the high-level clouds the 
greenhouse effect exceeds their albedo effect, which leads to a 
net heating in the lower atmosphere and an increase in the sur- 
face temperature (see also Fig . [7j» . 

To summarise, our parametrised cloud model is able to 
reproduce the mean Earth conditions, using measured cloud 
properties, cloud coverages, and the Earth global mean surface 
albedo. The Earth reference model can reproduce the mean Earth 
surface temperature and the Earth Bond albedo very well. The 
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Fig. 5. Radiative flux-pressure profiles of the Earth reference cal- 
culation including clouds. The upper diagram shows the down- 
ward radiative flux in the stellar (solid line) and the IR wave- 
length range (dotted line), the lower diagram the upward radia- 
tive flux. The position of the two cloud layers is denoted by the 
horizontal dashed lines. 



temperatures at the cloud positions also agree favourably with 
measured values. 

To mimic the climatic effects of clouds it is a common ap- 
proach to adjust the planetary s urface albedo in one -dimensional 
clear sky calculations (cf. e.g. lSegura et alj d2003l) ). While this 
approach can reproduce the correct surface temperature, these 
models have the shortcoming, amongst other things, that they 
are unable to reproduce the correct planetary Bond albedo. The 
parametrised cloud model, though, is able to reproduce both pa- 
rameters. 



5. Earth-like planetary atmospheres 

In this section we extend our studies of the climatic effects of 
clouds to Earth-like extrasolar planetary atmospheres at orbital 
distances for which the stellar energy input at the top of the at- 
mosphere equals the solar constant. The cloud parametrisation 
used is the same as for the Earth, but the cloud coverages of the 
high and low-level cloud layers are varied. To limit the number 
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Fig. 6. Surface temperatures of planets around the four different stars considered as a function of the cloud coverages of high and low- 
level clouds. Upper left diagram: F-type star, upper right diagram: G-type star, lower left diagram: K-type star, lower right diagram: 
M-type star. The solid lines on the contour surface and in the x-y-plane denote the physical limits of the cloud parametrisations. The 
parameters, for which a mean Earth surface temperature (288.4 K) is obtained are also marked by solid lines (see text for details). 
Each contour square side represents a change of 2.5% in cloud coverage. The dashed line in the x-y-plane separates the parameter 
regions in which partial overlap of the two cloud layers occurs or not. 



of free parameters the cloud overlap is calculated from the cloud 
coverages. As long as the sum of the coverages of both cloud 
layers is < 100%, no overlap occurs. If the sum of the two cloud 
coverages exceeds 100%, the cloud layers are partially overlap- 
ping with the smallest amount of overlap possible. As an exam- 
ple, a low-level cloud cover of 20% and high-level cloud cover 
of 80% results in no overlap, while for 60% and 80% cover- 
ages, the clouds are overlapping by 40%. Thus the resulting two 
dimensional parameter space, spanned by the two cloud cover- 
ages, is divided into two regimes of partially overlapping clouds 
and clouds without overlap (see Figs. l6l and [TOl. 

In order to quantify the climatic effects of clouds in Earth- 
like planetary atmospheres calculations have been carried out 
over the whole range of the two-dimensional parameter space 



of possible cloud cover combinations. For each stellar type more 
than 150 atmospheric models have been evaluated. The resulting 
surface temperatures and Bond albedos are presented in the next 
subsections. First implications of the position of the habitable 
zones of the four different central stars considered are discussed 
at the end of this section. 

5.1. Surface temperatures 

Figure [6] shows the surface temperatures as a function of the 
cloud coverages for the four different central stars. The various 
lines shown on the contour surfaces indicate different regions 
of physical conditions. The lowest line marks the validity range 
for the low-level cloud. For parameters below the correspond- 
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Fig. 7. Basic effect of each cloud type on the surface temperatures for the four different central stars. The upper diagrams show the 
surface temperature as a function of the cloud coverage of the high-level cloud (left diagram) and low-level cloud (right diagram). 
The lower diagram shows the derivative of the surface temperature with respect to the cloud coverage, i.e. the change of surface 
temperature with cloud coverage. 



ing line in the x-y-plane, the temperature of the low-level cloud 
drops below 260 K, which is the lower temperature limit for the 
approach (see Sect. [2]), i.e. it represents the freezing limit for 
water droplet clouds. The uppermost line on the other hand indi- 
cates the validity range of the high-level clouds. Above that line 
the high-level cloud temperature is higher than 260 K, where the 
ice cloud would liquefy. The middle line denotes the parameters 
for which a mean Earth surface temperature of 288.4 K occurs. 
It should be noted that all three lines become important only for 
the K and G-type stars. Due to the low resulting temperatures for 
the F-type star, the limit for the high-level cloud is not reached, 
while for the M-type star the temperature of the low-level cloud 
always stays above 260 K. 

Using the measured Earth surface albedo, results imply that 
none of the four clear sky model atmospheres achieves surface 
temperatures equal to the Earth mean value of 288.4 K0. Even 
if the incident stellar energy is the same for all types of stars and 
equals the solar constant, the resulting clear sky surface tem- 
peratures are obviously different (see Figs. [6] |7), which makes 
the influence of the wavelength dependence of the incident ra- 
diation quite clear, although no clouds are present and only the 
wavelength dependent opacities of the gas are causing the dif- 
ference in these cases. The surface temperature of a cloudless 
Earth-like planet orbiting the F-type star, for example, is lower 
than the measured mean Earth surface temperature, while the 
other surface temperatures are higher than 288.4 K. This is espe- 
cially the case for the M-type star, which has substantially more 
radiation flux in the visible and IR wavelength regions than the 
other three kinds of stars (cf. Fig. [3]). Consequently clouds are 
required to reach 288.4 K in all investigated situations. Whereas 
the cloud layers have to produce a net heating effect in the case 
of the F-type star, the albedo effect of clouds has to be stronger 
than the greenhouse effect for the other types of stars if mean 



1 This is of course the reason why the surface albedo has to be ad- 
justed to mimic the effect of clouds in clear sky model calculations for 
Earth. 



Earth surface temperatures are to be achieved under such condi- 
tions. 

In order to illustrate the basic climatic effect, the results of 
single cloud layer calculations are shown in Fig. [7] which sum- 
marises the corresponding slices through the 3D-temperature di- 
agrams of Fig. [6] As a rule the surface temperature decreases 
with increasing low-level cloud coverage, i.e. water droplet 
clouds exhibit a net albedo effect for all stellar spectra. The max- 
imum temperature decrease of about 30 K at full cloud cover is 
about the same for all types of central stars. The change of the 
surface temperature with increasing low-level cloud coverage is, 
however, not uniform. The most significant difference appears 
for the M-type star, while the for the F-type star there is almost 
no change in the temperature (see lower right panel of Fig. |7). 
This effect is caused by different absorption properties of the 
clouds at different wavelengths in conjunction with the incident 
stellar spectra. The M-star spectrum (Fig. [3} is more extended 
into the IR where the frequency dependent spectral energy can 
be partly absorbed by the low-level water cloud (cf. Fig. [TJ. In 
contrast to the M-type star, the spectrum of the F-type star has its 
maximum flux at a wavelength range around 0.5 /an, where the 
water cloud has an almost constant scattering albedo of about 1 
(see Fig. [T}, which means that scattering dominates the radiative 
transfer, and the remaining absorption has only a minor influence 
on the temperature. 

Due to their net greenhouse effect the high-level ice clouds 
cause an increase in the surface temperature with increasing 
cloud cover in all cases (see Fig. |7). The maximum of the tem- 
perature increase caused by the greenhouse effect, however, de- 
pends on the properties of the central star considered. The albedo 
effect depends directly on the incident stellar spectrum, but the 
greenhouse effect is determined by the thermal emission of the 
lower atmosphere, which itself is an indirect consequence of the 
stellar radiation. For example an Earth-like planet around the M- 
type star has a maximum greenhouse effect of about 1 1 K at 
100% high-level cloud cover compared to the clear sky case. 
In the case of the F-type star the greenhouse effect results in 
an increase of surface temperature causing a temperature alter- 
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Fig. 8. Illustration of the efficiency of the greenhouse effect 
at different atmospheric temperatures of the layer below the 
high-level cloud. The solid line represents the absorption opti- 
cal depth of the high-level ice cloud. The black-body radiation 
fluxes for planets around the different central stars (F-type star: 
dashed-dotted, G-type star: dotted, K-type star: solid, M-type 
star: dashed) were scaled for comparison. The temperatures were 
calculated with 50% high-level cloud cover. The vertical dotted 
lines denotes the positions of the maxima of the black-body ra- 
diation fluxes. 

ation of about 20 K. In principle the greenhouse effect becomes 
smaller for larger temperatures in the lower atmosphere (Fig.|7j 
left panel). The effectiveness of the greenhouse effect depends 
directly on the absorption characteristics of the high-level cloud 
in the thermal wavelength range in conjunction with the temper- 
ature dependent infrared emissions of the lower atmosphere. 

To illustrate the absorption characteristics of the clouds as a 
function of different atmospheric temperatures, we assume that 
the transmission of the atmospheric layer below the cloud layer 
can be represented by black-body radiation with the correspond- 
ing atmospheric temperatures. Figure[8]shows the absorption op- 
tical depth of a high-level ice cloud with 50% coverage in com- 
parison to black-body radiation fluxes of different atmospheric 
temperatures derived from model calculations for the four differ- 
ent central stars. Compared with the other Earth-like cases, the 
planet around the F-type star has the lowest atmospheric temper- 
atures, which results in the biggest shift of the black-body radi- 
ation to longer wavelengths. Since the absorption at the position 
of the maximum of the black-body radiation is more than 12% 
smaller for the M-star than for the F-star, the greenhouse effect 
is much more significant in the latter case. Thus one obtains a 
temperature-induced wavelength shift for the absorbable radia- 
tion, which results in increasing absorption and a more effective 
greenhouse effect. As the temperature in the lower atmosphere 
increases with larger high-level cloud cover, the maximum of 
the black-body radiation is more and more shifted to smaller 
wavelengths which in turn decreases the efficiency of the green- 
house effect. The amount of the greenhouse effect is in that way 
self-limited. This can clearly be seen in Fig.|7](left lower panel) 
which shows the flattening of the temperature gradients at high 
ice cloud coverages. 

The climatic influence of clouds in Earth-like planetary at- 
mospheres becomes much more complex for a multi-layered sit- 
uation, where the greenhouse effect and the albedo effect of the 
high and low-level clouds are interacting non-linearly, resulting 
in the 3D-temperature planes displayed in Fig. [6] The efficiency 
of the ice cloud greenhouse effect, for example, is much more en- 
hanced if a low-level cloud is present below the high-level cloud. 
The albedo effect of the low-level cloud is compensated by the 



310 



_ 




• • 


--• " . • 
- -' ' 


F-tvnp ^tar 




G-type star _ 




K-type star 


= 1 1 1 h 


M-type star 




' 1 1 1 1 = 

■ ■ ■ F-type star - 




G-type star . 


- 


,. K-type star 




_ .. - " "\ -M-type star 







20 40 60 SO 100 

high-level cloud cover x (%) 

Fig. 9. Effects of multi-layered clouds on the surface temper- 
atures for the four different central stars. The upper diagram 
shows the surface temperature in dependence of the high-level 
cloud coverage with 100% low-level clouds. The lower diagram 
shows the derivate of the surface temperature with respect to the 
cloud coverage. 

ice cloud greenhouse effect. Therefore the resulting temperature 
difference for increasing high-level cloud coverage is much more 
pronounced at 100% low-level cloud cover than for a single ice 
cloud layer. This can clearly be inferred from Fig. [6] for all four 
considered scenarios. 

The enhancement of the greenhouse effect is, of course, a 
result of lower atmospheric temperatures, reduced due to the 
presence of the low-level cloud and its albedo effect. Due to the 
reduced temperatures the intensity maximum of the absorbable 
radiation is shifted to longer wavelengths, thereby increasing the 
efficiency of the greenhouse effect (cf. Fig. [8). The increase of 
the surface temperature at full low-level cloud cover with in- 
creasing high-level cloud coverage is shown in detail in Fig. 
[9] which combines four different slices through the 3D-surface 
temperature diagrams of Fig. |6l respectively. Obviously, the 
temperature increase becomes weaker at large ice cloud cover- 
ages for all cases, as also revealed by the maxima of the cor- 
responding temperature gradients (F-type star: as 60%, G-type 
star * 50%, K-type star * 40%, M-type star ~ 40%) shown in 
the lower panel of Fig. [9] This is caused by the change of the 
energy transport mechanism from radiative transfer (at smaller 
high-level cloud coverages) to convection (for larger high-level 
cover) at the calculated heights of the upper cloud layers in the 
planetary atmospheres. The switch from convection to radiative 
transfer is also responsible for the second feature in the tem- 
perature slope of the M-star planet occurring at low ice cloud 
coverages. However, note that the highest (lowest) values of 
the surface temperatures affected by clouds are reached for sin- 
gle ice clouds (low-level clouds) with 100% coverage (Fig.[6])Q 

2 The applicability limitations of the cloud description should thereby 
be kept in mind. 
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Fig. 10. Planetary Bond albedo of planets around the four different stars as a function of the cloud coverages of high and low-level 
clouds. Upper left diagram: F-type star, upper right diagram: G-type star, lower left diagram K-type star, lower right diagram M-type 
star. Each contour square size represents a change of 2.5% in the cloud coverage. The dashed line in the x-y-plane separates the 
parameter regions in which partial overlap of the two cloud layers occurs or not. 



The difference between the maximum and the minimum sur- 
face temperature again depends on the characteristics of the cen- 
tral star. The F-type star causes the largest temperature variation 
(x 50 K), whereas the temperature difference for the M-type star 
is the smallest (« 40 K). 

We note that the mean Earth surface temperature, which re- 
quires the presence of clouds, is not uniquely obtained by one 
single combination of low and high-level cloud coverages. On 
the contrary, several cloud cover combinations can in principle 
result in a surface temperature of 288.4 K as indicated by the cor- 
responding contour lines in Fig. [6] However, the range of cover- 
age combinations for which mean Earth surface conditions can 
be reached differs between the stellar types from the M-type star 
(smallest parameter range) to the F-type star (largest parameter 
range). Due to its higher incident stellar flux in the visible and IR 
wavelength region, the M-type star requires a high cover (about 



60%) of low-level clouds as a minimum to achieve a cooling to 
288.4 K. For the F-type star on the other hand only about 10% of 
high-level clouds are needed for the required additional heating. 

5.2. Bond albedo 

The Bond albedo is generally an important quantity to charac- 
terise the influence and properties of clouds in planetary atmo- 
spheres. Here the Bond albedo values have been obtained for the 
whole two-dimensional parameter space of cloud coverages. In 
Fig.[l0]the calculated planetary Bond albedos are shown for all 
four central stars considered. As expected, the albedo increases 
regardless of the cloud type with increasing cloud cover; i.e. in 
absence of any cloud the Bond albedos are given by the proper- 
ties of the gas only and have very low minimal values between 
0.07 for the M-type star and 0.17 for the F-type star. Overall, 
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the albedo values for an Earth-like planet around an M-type star 
are lower than for the other stellar types, whereas for the F-type 
star the Bond albedos are the highest. This is the result of direct 
absorption of the incident stellar radiation by clouds and gas. 
Since the clouds can absorb a larger fraction of the M-star spec- 
trum (see Figs. Q] and [3), less of that light is reflected back into 
space, resulting in a lower Bond albedo. 

The low-level clouds have a larger impact upon the Bond 
albedos than the ice clouds because their optical depth is twice 
as large in the wavelength range of the stellar spectra (cf. Figs.Q] 
and 0. Therefore the Bond albedos caused by low-level clouds 
alone are higher than the albedos of single ice cloud layers. One 
would expect that the Bond albedos become maximal in the case 
of total cloud cover of both cloud layers. This is the case for 
the F-, G- and K-type conditions, whereas for the M-type star 
the single low-level cloud layer with 100% coverage yields the 
highest Bond albedo. This is related to a slight anomaly which 
is evident for the M-type star and also (but less pronounced) for 
the K-type star (Fig.lTOt. At total low-level cloud cover the Bond 
albedo starts to decrease with increasing high-level cloud cover- 
age. Where the ice cloud cover exceeds about 50% the albedo 
values increase again. 

This effect can be illustrated by the radiative flux-pressure 
profiles for different high-level cloud coverages and full low- 
level cloud cover in the atmosphere of an Earth-like planet or- 
biting an M-star (Fig. ITTb . The ice cloud presence affects the 
downward radiative flux as expected (upper panel of Fig. fTTT i by 
enhanced absorption and scattering capability. But above the up- 
per cloud level the upward radiation flux for full ice cloud cov- 
erage is higher than the corresponding radiation flux for only 
50% cover (see lower diagram of Fig.lTO. This can be explained 
by increased gas absorption of the incident stellar radiation for 
these specific atmospheric conditions. This effect deserves fur- 
ther investigation especially in view of improved gas opacities. 
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Fig. 11. Radiative flux-pressure profiles of different model calcu- 
lations for the M-type star to illustrate the albedo anomaly. The 
calculations have been performed with a 100% cover of low- 
level clouds and three different high-level cloud coverages. The 
upper diagram shows the downward radiative flux in the stellar 
wavelength range, the lower diagram the upward radiative flux. 
The position of the two cloud layers is denoted by the horizontal 
dashed lines. 



5.3. Habitable zone positions influenced by H 2 clouds - first 
implications 

The potential habitability of terrestrial planets depends on 
their surface conditions, especially on the surface temperatures. 
Usually the possible existence of liquid water on the planetary 
surface is considered as an indication for habitable conditions. 
Below we refer to the measured mean Earth surface temperature 
of 288.4 K as characteristic for the positions of habitable zones 
around different kinds of central stars. 

As already pointed out in Sect. l3.3l the orbital distance of the 
Earth-like planets to their host stars have been determined here 
in a way that the incident stellar energy matches the solar con- 
stant regardless of the atmospheric details (see Table [3]). Using 
these distances, the calculated planetary surface temperatures do 
not comply with the mean Earth surface temperature assuming 
clear sky conditions as discussed in detail in Sect. 15. II Habitable 
conditions can yet be achieved even for these positions due to the 
complex climatic effects caused by H2O cloud layers as for the 
Earth for example. As previously described, mean Earth surface 
temperature conditions can in principle be reached for several 
cloud coverage combinations. Adjusting the stellar energy input 
at the top of the atmosphere by changing the distance between 
the central star and planet properly can, however also result in 
Earth-like conditions for clear sky atmospheres. The deviations 
between these approaches can clearly be seen in Fig. [12] 

Using clear sky atmospheres the orbital distances of the plan- 
ets have to be increased (by about 2% for the M-type star and 
3% in the case of G-type and K-type stars) to compensate for 



the missing net cooling effect of the H2O cloud layers and to 
achieve Earth-like conditions on the surface. For the F-type star, 
the distance needs to be decreased by about 3%, because in 
the clear sky case the surface temperature is below 288.4 K. 
Consequently, planets with a measured mean Earth cloud cover 
can be located closer to the central star than planets with clear 
sky atmospheres with the exception of the F-type star. 

The p ositions of Earth-like plane ts in the habitable zones cal- 
culated bv lSegura etafl d2003l |2005|) are also indicated in Fig. [HI 
for comparison. These distances have been derived from clear 
sky calculations for the Earth atmosphere with an adjusted plan- 
etary surface albedo of 0.2 to account for the already mentioned 
net cooling effe ct of cloud and is used f or all different types 
of central stars. ISegura et all d2003l [2005) had additionally to 
change the distances of these planets to their host stars to reach 
a temperature of 288 K at the planetary surface even with the 
modified surface albedo. These positions deviate from the dis- 
tances of planets with an Earth-like cloud cover (this work) be- 
tween 17% in case of the K-type star, about 10% for the F-type 
star and less than 1 % for the M-type star. These differences are 
partly the consequence o f the increased surface albedo used by 
ISegura et all d2003l 120051) . which is 54% larger than the mea- 
sured Earth mean value used here. A larger surface albedo leads 
to less absorbed stellar flux at the surface. Therefore the planets 



3 This tuning, however, is only possible if the target surface temper- 
ature is prescribed, while in case of the model including clouds the sur- 
face temperature is a result of the calculations. 
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Fig. 12. Positions of the habitable zone around different types 
of central stars affected by clouds. Squares mark the distances at 
which the incident stellar flux matches the solar constant. Circles 
indicate the distances for which a mean Earth surface tempera- 
ture is achieved using the measured Earth cloud cover. The tri- 
angles show the positions for planets with clear sky atmospheres 
and a mean Earth surfa ce temperature. The corresponding dis- 
tances derived by Segu ra et alj ([2003, 2005) using a clear sky 
model with adjusted surface albedo are marked by diamonds. 
The maximum effect of the clouds on the distances are marked 
by arrows. 



in the work of Segu ra et~a"fl ([2003, 2005) are located closer to the 
central star (except for the M-type star) compared to our findings 
including the effect of clouds. Still, it should be noted that other 
differences between the two models do also affect the deviation 
in the positions of the habitable zones. An important difference 
in the modelling approaches is certainly the different treatment 
of the atmospheric chemistry. Here Earth-like chemical p rofiles 
are used (see Fig. [2]), whereas Se gura et"al] d2003l 120 05) deter- 
mined the chemical composition from a photochemical model. 
The obvious deviation in the case of the K-type star is mainly 
caused by the different stellar input spectra used (cf. Sect. 13.31 ). 
But apart from the well known overall effect that the positions 
of the habitable zone gets closer and clos er to the central star b y 
changing its spectral type from F to M (cf. Kasti ng et alj (Q993)), 
the special spectral characteristics of the stars are affecting the 
positions of the habitable zones differently due to their complex 
interaction with the atmospheric cloud layers (see Fig. [121. 

The maximum effect of mean Earth clouds on the position of 
the habitable zone are denoted by arrows in Fig. Q~2] These dis- 
tances have been derived with 100% of either cloud type while 
still obtaining a mean Earth surface temperature of 288.4 K. 
Using low-level clouds for a maximum cooling effect, the plan- 
ets can be located up to 15% closer to the central star com- 
pared to a clear sky planet. The distances can be increased by 
up to 35% when sin gle high-level cloud l ayers are present. The 
distances derived by Segura et al. ( 20031 120051) are in all cases 
within these boundaries, except for the already discussed K-type 
star. 

It is important to note that the distances derived in this work 
are only indications of the position of the habitable zone, they 



do not indicate its boundariefl The inner boundary is deter- 
mined by the runaway greenho use effect while the outer edge is 
determ ined by CO2 clouds (cf . Kastingl d!988l) and ISelsis et all 
(120071) ). Since these kinds of atmospheres are in general far from 
being Earth-like, our parametrised cloud description, which is 
based on measurements in the Earth atmosphere, is unsuitable 
to study the effects of clouds on the boundaries of the habitable 
zone. Further studies of the effects of clouds on habitable zones 
require therefore a cloud microphysics model to evaluate the 
cloud properties (e.g. optical thickness, wavelength dependent 
optical properties, altitude) properly, since measurements and 
therefore also parametrised descriptions of them are not avail- 
able for these more extreme physical conditions. A more detailed 
analysis of the influence of clouds on the position and extension 
of the habitable zone will be addressed in a forthcoming publi- 
cation. 



6. Summary 

We developed a parametrised cloud description scheme for 
Earth-like planetary atmospheres based on measurements of 
clouds in the Earth's atmosphere. The optical properties of the 
clouds were calculated with the Mie theory combined with 
an equivalent sphere approach for the ice crystals. This cloud 
scheme was coupled with a one-dimensional radiative convec- 
tive climate model, and its applicability was tested on an Earth 
reference model. Model calculations were made for Earth-like 
planetary atmospheres of planets orbiting different types of main 
sequence dwarf stars: F, G, K, and M-type stars. 

It was shown that the albedo effect is only weakly depen- 
dent on the incident stellar spectra because the optical proper- 
ties (especially the scattering albedo) remain almost constant 
in the wavelength range of the maximum of the incident stel- 
lar radiation. The greenhouse effect of the high-level cloud on 
the other hand depends on the temperatures of the lower atmo- 
sphere, which in turn are an indirect consequence of the differ- 
ent types of central stars. The efficiency of the greenhouse ef- 
fect increases with smaller temperatures in the lower atmosphere 
and decreases with higher atmospheric temperatures. As a rule 
the planetary Bond albedos increase with cloud cover of either 
cloud type. However, for the K and M-type star an anomaly 
was found resulting in a decreasing Bond albedo with increas- 
ing cloud cover for a certain region in the parameter space. This 
anomaly is caused by enhanced gas absorption under these spe- 
cific atmospheric conditions. 

Planets with Earth-like clouds in their atmospheres can be 
located closer to the central star or farther away compared to 
planets with clear sky atmospheres. The change in distance de- 
pends on the type of cloud. In general, low-level clouds result in 
a decrease of distance because of their albedo effect, while the 
high-level clouds lead to an increase in distance. 

Apart from these climatic effects, clouds also affect the plan- 
etary emission and transmission spectra, which is important con- 
cerning the detectability of spectral signatures of e.g. biomarker 
molecules. The investigation of the influence of clouds on these 
effects will be done in further studies. 
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4 See e.g. ISelsis et all (12007) for a study on the boundaries of the 
habitable zone of planets around different types of central stars. 
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